A fractal approach to the rheology of concentrated cell suspensions 
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Results on the rheological behavior of novel CHO cell suspensions in a large range of concen- 
trations are reported. The concentration-dependent yield stress and elastic plateau modulus are 
formalized in the context of fractal aggregates under shear, and quite different exponents are found 
as compared to the case of red blood cell suspensions. This is explained in terms of intrinsic micro- 
scopic parameters such as the cell-cell adhesion energy and cell elasticity but also the cell's individual 
dynamic properties, found to correlate well with viscoelastic data at large concentrations (0 > 0.5). 



The rheology of complex fluids has been studied ex- 
tensively over the past decades [1] and has revealed very 
intriguing behaviors, in particular properties of suspen- 
sions, either micronic or colloidal, are still a subject of 
interest [H, 0, 0, @. Classical behaviors of suspensions 
reveal shear-thinning effects usually, but other unusual 
ones like shear-thickening Q (i.e. viscosity increase with 
shear rate) or yield stress have been observed [2|, l5[ . The 
yield stress is the critical value of the shear stress needed 
to induce flow for a given fluid. It is closely related to the 
internal structure of the fluid therefore its ability to form 
(or break) particle clusters under flow. In this respect 
many studies have focused on solid sphere suspensions. 

On the other hand, there are much less works dedicated 
to suspensions of deformable particles, such as biological 
cell suspensions. The main works can be found in the 
field of blood rheology. Suspensions of Red Blood Cells 
(RBC) within plasma were first investigated by Chien 

0. 01 and revealed a shear-thinning behavior, but a more 
detailed inspection of the viscosity-shear rate diagrams 
showed that at low shear rates, the stress level is close to a 
constant <j s (Pa), called the yield stress. The well-known 
Casson's model |9[, y/a = \f& s + y/JPy, relating the shear 
stress a to the shear rate 7 (/i being a constant viscosity) 
can be used to determine the yield stress. Chien and co- 
authors obtained <j s for a large range of hematocrit (H), 

1. e. the RBC volume concentration [TO]. They showed a 
relationship of the type a s ~ (H — b) 3 (b being a constant 
hematocrit) as also observed in a recent work [111 ]. 

It is still not known yet whether this type of behavior 
is universal, or rather it could depend on cell type, cell 
shape or other biological effects such as cell adhesion or 
cell elasticity. In particular, one proposed explanation of 
the yield stress in RBCs suspensions is based on the ex- 
istence of " rouleaux" which build due to cell interactions 
and exhibit large shape aspect ratios [8] and a fractal 
dimension D. Therefore it is necessary to apply strong 
enough stresses in order to break such aggregates, in close 
relation with the yield stress. 

In this letter we propose to investigate the rheology of 
a new cell suspension, consisting of CHO cells (Chinese 
Hamster Ovary cells) in a large range of concentrations. 
Such cells are commonly used in biology, easy to culture, 



and can be genetically modified to induce different adhe- 
sive properties. These cells are spherical when suspended 
in a culture medium, and organized in a specific manner 
leading to particular aggregation patterns of fractal type. 
This leads to the determination of scaling laws based 
on fractal exponents (for the yield stress <j s and elas- 
tic modulus Go) which are seen to be non universal but 
dependent on cell type. The flow curves constitute a ba- 
sis to test classical empirical models (Bingham, Casson, 
Herschel-Bulkley models) and other ones [H, El based 
on kinetic theories describing the rupture and formation 
of particle clusters. The latter ones successfully relate 
macroscopic effects to microscopic parameters, such as 
the cell-cell adhesion energy and the cell elasticity. These 
microscopic parameters remarkedly match those found in 
the litter ature using other techniques. This is important 
in the context of recent studies related to tumour growth 
EH EH, EH which consider cell assemblies with interac- 
tions as well as cell elastic deformations. Furthermore, 
this study emphasizes the relationship between the dy- 
namic rheological properties of suspensions [17[ and the 
single cell properties. 

Experiments - In our model system, CHO cells are 
grown in culture medium using standard T75 boxes un- 
der proper conditions (37°C, 5%C02), until they are 
at 70% confluence, then are detached using trypsin, are 
mixed together and centrifugated at 1200 rpm, a high 
enough velocity to get a concentrated suspension, but 
slow enough in order to maintain the cells alive. Cell vol- 
ume concentration <j) (i.e. similar to the hematocrit H) 
is determined accurately after centrifugation in hemat- 
ocrit tubes containing the CHO cells. Then the right 
amount of remaining supernatant is removed until the 
desired concentration is obtained (between and 60%). 
Different experiments were carried out on a conventional 
Rheometer (Bohlin Gemini 150). Both steady shear and 
oscillatory measurements were made at T = 20°C. Due 
to the large amount of cells needed (we usually require 
twelve T75 flasks in order to obtain a volume of roughly 
0.3 mL of cells), we chose to use a plate-plate geome- 
try (20 mm diameter) with a small gap (between 400 [im 
and 1 mm) for the concentrated suspensions whereas the 
smaller concentrations (below 10 %) were tested using the 
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FIG. 1: Viscosity 77 (Pa.s) vs. shear rate 7 (s 1 ) at different 
volume concentrations 6 from to 60% 
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FIG. 2: Shear moduli G' and G" (Pa) vs. angular frequency 
cj at different volume concentrations ranging from 20% to 
60% 



60mm cone-plane geometry (2° angle). Typically in our 
fluid, the suspended cells are spherical and monodisperse 
with a radius a ~ 10 fim. 

Results - Experimental results for constant steady 
state shear rate 7 are presented in FigHJ The viscosity 77 
is shown to vary over several decades, within shear rates 
typically between 10 _3 s _1 and 10 3 s _1 . In some cases, we 
limited ourselves to the higher shear rates because of ex- 
perimental reasons (i.e. steady state not reached). By a 
first inspection of the curves, we recognize the signature 
of a yield stress fluid as depicted by the slope close to 
— 1 in the viscosity-shear rate diagram (or equivalently a 
constant shear stress at low shear rates), especially at the 
largest concentrations 0, which will be particularly of in- 
terest here. The existence of this yield stress is attributed 
to weak interactions which can exist after preparation of 
the system. Already existing proteins are available on cell 
membranes and can be recruited to form bonds, leading 
to particular structure arrangments. This explains the 
presence of a yield stress related to the formation of such 
structures. The yield stress is found to depend on volume 
concentration in a manner to be discussed later. 

A second series of experiments was carried out in order 
to study the systems under oscillatory strains at angular 
frequency uo. Small deformations (1% or less) within the 
linear regime were performed in order to characterize the 
elastic modulus G'(uj) and the loss modulus G"(uj). We 
find an interesting behavior as shown in Fig|2j Moduli 
G" (co) usually prevails over G' (uo) at small concentrations 
(e.g. (j) = 0.2), but as <p increases, the system becomes 
elastic with a much larger G'(uj). This behavior is the 
signature of a viscoelastic medium, due to the fact that 
interactions between elastic cells become effective at large 
concentrations ((j) > 0.4). The slow increase of the elastic 
modulus G' against frequency reveals the presence of a 
so-called 'elastic plateau' modulus (Go) determined by 
the value of G'(cu) at intermediate rates (lrad/s typi- 
cally). The presence of elasticity has been observed pre- 
viously for RBC suspensions [18], above a critical vol- 



ume fraction around <p = 0.2, and is believed to come 
from the elasticity of the cells as they are packed more 
closely at large concentrations such as the ones also en- 
countered in tumour spheroids [l5[. Finally, we observe 
that the trends in the G' — G" plots for large concen- 
trations (<j) > 0.5) are remarkedly similar to previous mi- 
crorheological results obtained on single cells pjl [10, Hl| . 

As in the case of suspensions, we define a maximum 
packing fraction </>o (which is usually 0.64 or even 0.74 for 
solid spheres in a face-centered-cubic crystal) , depending 
on cell elasticity, i.e. their compactness [22]. Due to 
the presence of soft spherical cells, it is expected that 
the value of 0o will be in this range. </>o is determined 
using the reduced viscosity plot ^ as a function of <fi (77 
at a shear rate of 10 3 s -1 , 770 = 0.0014Pa.s the solvent 
viscosity). In our case, this data (not shown) is found 
to match the well-known equation ^- = (1 — A)- 2 -5</>o 
proposed earlier [23[, this providing the value </>o ~ 0.65. 

In order to investigate the effect of the volume con- 
centration (/>, we first need to obtain the flow curve (7(7) 
of the suspensions, as well as the relevant parameters, 
such as the yield stress a s . From the viscosity curve in 
FigHJ we plot the stress a = 77(7)7 vs - shear rate 7 
in Figl3l and fit the data with the Herschel-Bulkley law 
[5[ a = a s + i^7 n , where K is a constant, and n is a 
shear-thinning exponent ranging between and 1 (1 is 
for a Bingham fluid, and the case of the Newtonian fluid 
is recovered for n = 1, cr s = 0). Parameters have been 
optimized using a standard Newton-Raphson method. 

This leads to the determination of the yield stress cr s 
as a function of volume fraction (j). Such measurements 
are usually difficult [10[ because of possible slip, sedi- 
mentation and evaporation [11, 24]. Care has been taken 
to avoid such problems, therefore only shear rates larger 
than 10 _3 5 _1 (lowest value) are considered. The empir- 
ical Herschel-Bulkley model (involving a yield stress) is 
then used when sufficient data points are available. The 
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FIG. 3: Determination of the yield stress <j s using Herschel- 
Bulkley's model. 




FIG. 4: Yield stress a s and Shear elastic modulus Go = 
G'(u) — lrad/s), vs. volume concentration 0, log-scale. 



fits are in satisfactory agreement with the data which 
gives good confidence in the values of the yield stresses 
for cf) > 0.42. Another attempt has been made using Cas- 
son's model and gives similar data. The Bingham model 
seems to give less accurate values. 

The values of the yield stresses as a function of volume 
concentration are plotted in FigH^, to be compared with 
the evolution of the shear plateau modulus Go (value 
of G' at a typical angular frequency uo = Irad/s) in 
Figj4)3. These two plots show power law dependences of 
the form o~ s ~ (j) mi and Go ~ </> m2 and remind previous 
results [12[ obtained in the case of the rheology of RBCs 
suspensions, at least for the yield stress <j s . From FigH^- 
b we find that mi ~ 8.4 and vrt2 ~ 11.6. These exponents 
are quite different from the ones obtained in the case of 
RBCs suspensions and this will be discussed below. 

Modelling - As seen above, rheological modelling of 
such suspensions should therefore predict shear-thinning 
behavior, as well as yield stress properties at low shear 
rates 7^0 and a concentration dependence of <j s and 
Go- In addition, cell suspensions correspond to aggre- 
gated systems (see FigJSJ). Under flow, their structure 
is based on the persistent remodelling of the cells with 
respect to each other as they exhibit deformations, ro- 
tations, possible rolling and/or separation. During such 
events, cells may form clusters of size Rf to be com- 




FIG. 5: Phase contrast microscopy of CHO cell suspension : 
10% and 52%. Same scale for both images. 



pared with the cell size a (radius). The formation and 
destruction of cell clusters is the major ingredient to un- 
derstand the rheological properties of the cell system, in 
order to explain our previous data. For example, when 
sheared under stress <r, clusters break into smaller ones, 
leading to shear-thinning effects. On the other hand, the 
possible encounter of clusters leads to the formation of 
larger structures, increasing the viscosity. Using a se- 
ries of cliches like in Fig El we measure the size Rf of 
clusters at rest, found to vary as a power-law of the 
type ^ ~ N 1 / [25], where N is the number of cells 
in a cluster and D the fractal dimension. For our case, 
we find that D ~ 1.47 for 2D images. Thus, in three 
dimensions, we expect a fractal dimension of the order 
D ~ 2 (26[. This number is similar to the ones found for 
RBCs suspensions, although the scaling exponents for 
yield stresses are quite different. 

In the semi-empirical model proposed by Snabre and 
Mills [H, Qjl, the formation and dissociation of clus- 
ters under flow is taken into account. A change in Rf 
as a function of the applied shear stress is assumed : 
= 1 + (^-) m , where m is a dimensionless parameter, 
a* is a critical stress related to the interfacial adhesion 
between cells : a* = T/a, and T is the cell adhesion free 
energy. Using the concept of effective medium with vol- 
ume fraction 4>a = </> (^-) 3-jD , one assumes an effective 
viscosity n (a) = rjo 1 ~j^ 4 , and obtains the constitu- 

tive equation [27| which contains the yield stress given 
by a s ~ <t*(^-) ™(3--o) m xhe last parameter to be used in 
the formula, (/>o, is the maximum packing concentration 
found previously. 

We use the previous model to explain our experimental 
data. The exponent mi =8.4 found for the yield stress <j s 
is plugged into the previous scaling law for determination 
of the parameter m = 0.078. This is smaller than the val- 
ues of m found for RBC suspensions (typically m ~ 0.3). 
This means that the size of clusters is not so sensitive to 
the applied stress, indeed one can consider that the cell 
aggregates are easier to form because of the round shape 
of the cells, in contrast with RBCs which need to bind 
in a very special way to form 'rouleaux', thus the stress 
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has more effect on the latter ones. We have obtained the 
value of the critical stress <r* = lAN/m 2 , and a corre- 
sponding value of T = 1.410 _5 7V/m. This value of a* 
is higher than the ones obtained for RBCs [l2[ but the 
interfacial energy V is in the range of the small values 
indicated for vesicles [28]. This is in favor of the initial 
assumption that few adhesion molecules are involved in 
the region of contact between the cells. 

Finally, we postulate a similar relationship [25[ for the 

shear elastic modulus Go ~ G*(^-) ™(3-°) , where G* is 
an effective elastic modulus, but we include an addi- 
tional exponent n to be determined. We come up with 
n = 0.056 and G* = 234Pa. This value of the reference 
modulus G*, as explained in the concept of fractal ex- 
ponents [25[, is to be related to typical values for single 
cells. In particular, it gives a Young's elastic modulus of 
E* = 702Pa (assuming that the cell is incompressible) 
which is typical for adherent Wild Type CHO cells, of 
the order 0.5 - 1 kPa as measured by AFM [1|[30[. 

To sum up, the system studied here provides unique 
features important for the rheology of biological suspen- 
sions and tissues. These concentrated cell suspensions 
behave as yield stress fluids (also called visco-plastic ma- 
terials), for which a fractal approach has been used. Un- 
der shear, the fractal structure changes and can be mod- 
elled using a yield stress a s and elasticity modulus Go 
related to the fractal dimension D. Two other parame- 



ters of interest have been introduced in the model: the 
adhesion energy T found in the range of typical values for 
cell adhesion energies (r ~ 10 -5 N/m), and an effective 
elastic modulus (E* ~ 700 Pa) in agreement with previ- 
ous microrheology experiments. We also found a similar 
behavior between the dynamic shear moduli G'{uj) and 
G"{uj) in this study (at large (/>), and the ones obtained 
from microrheological studies on sing le cells [H, H| El 
using various techniques. Both show slowly increasing 
dynamic moduli in terms of frequency, with the same 
relative positions. This idea probably deserves more at- 
tention and should be tested in the future, in particular 
further works may focus on the characterization of other 
cellular suspensions including cells with different elastic 
properties. 

Finally, such a study can naturally lead to the under- 
standing of biological tissues, by including stronger ad- 
hesion properties between the cells, or by taking into ac- 
count the addition of Extra- Cellular Matrix components. 
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